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^^ We propose using the Skyrme model on a lattice as an effective field theory of 

, meson— baryon interactions. To this end we construct a local topological density 

^ ^ ■ that involves the volumes of tetrahedra in the target space S and we make use 

of Coxeter's formula for the Schlafli function to implement it. We calculate the 
mean-square radius of a skyrmion in the three-dimensional Skyrme model, and 
find some surprises. 
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1 Why a lattice Skyrme model? 

The Skyrme model is a theory of a scalar field U{x) G SU{2) with the action 



yI Tr \d,U\' - ^-^Tr {[d,UU\dMU^') 



(1) 



The action possesses an SU{2) x SU{2) chiral symmetry which is sponta- 
neously broken. The Goldstone bosons are taken to represent pions, massless 
k\ I unless we add a symmetry-breaking term to S. The theory also contains soli- 

<^ . tons, stabilized by the 4-derivative Skyrme term shown, that have the prop- 

erties of baryons. The model furnishes a rich phenomenology of pion-nucleon 
interactions at low momenta. 
1^^ ' The Skyrme model does not really exist as a continuum theory, since the 

S . action is non-renormalizable. For this reason, existing treatments of the model 

are semiclassical, quantizing only the collective degrees of freedom of the soli- 
ton. Full quantization of the theory requires a cutoff. In considering the sector 
without skyrmions, this cutoff can be removed order by order in perturbation 
theory at the price of an ever-lengthening list of higher-dimension counter- 
terms. In the skyrmion sector, however, even this is difficult to accomplish. 
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We suggest Id that the hmitations of the "continuum" model be turned to 
advantage by treating the Skyrme model as an effective field theory without 
removing the cutoff. Now the cutoff is part of the specification of the theory. 
The form of 5* may be chosen freely, aiibject only to phenomenological tests. 

To write a lattice Skyrme model,au we first choose an action. We keep 
the cutoff — the lattice spacing — fixed. We apply the full panoply of lattice 
methods to calculate quantities of interest, going beyond semiclassical meth- 
ods and beyond perturbation theory. Renormalization consists of adjusting 
the bare couplings and the lattice spacing to match physical parameters. No 
continuum limit is necessary — or possible. 

2 Lattice topology 

Continuum solitons are classified by the winding number of the field configu- 
ration that maps the compactified 3-space {R^ + cxi} to the field space SU{2), 
both of which arc topologically equivalent to the 3-sphere 5''^. The differential 
volume in the target space is given by 

dT - ^e,jkTr {U-^d,U)iU-'d,U)iU-'dkU)d\ (2) 

which integrates to an integer n. On a lattice, we calculate this differential 
volume directly as follows. 

1. We cut each lattice cube into five tetrahedra: 




2. Each tetrahedron maps to a tetrahedron in S^. Oddly enough, there 
is no simple formula for the volume of this curved tetrahedron.^] Our 
solution to this problem begins by dropping perpendiculars to cut the 



"^Cf. Girard's theorem for a triangle in 5^, which says that the area is equal to the angular 
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curved tetrahedron into six quadrirectangular tetrahedra (qrt's), of which 
we show two: 




In general, a qrt can be constructed by drawing three mutually perpen- 
dicular line segments P0P1P2P3 and connecting all the vertices: 




(This works in S^ as well as well as in R^.) All the faces of a qrt are 
right triangles. Three of the dihedral angles of a qrt are right angles; the 
other three may be labeled a, (3, 7. 

Polytopcs in S*" were studied by Schlafli in the 19th century, and in 1858 
he wrote down a differential equation for the volume of a qrt as a function 
of a,P,"f. This equation was solved by Coxctcr in 1935. The volume V 
is given by the Schlafli function S according to 



1 /TT 7T 

V{a, (i, 7) = -S* (^- - a, /3, - - 7 



S{x,y,z) = ^ 



D — sin X sin z \ cos 2mx — cos 2toj/ 
D + sin X sin z 



(3) 

cos 2mz — 1 



(4) 
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where D = v sin^ a s\v? 7 — cos^ j3 is the "angular excess" — which van- 
ishes for a EucHdean qrt. In fact y ^ as I? ^ 0, as expected. 

3 The lattice action 

The topological density defined above will not be meaningful until we resolve 
ambiguities associated with what is known as "topology dropping through the 
lattice." A field configuration U{x) specifies the vertices of all the tetrahedra 
in S^ . Given the vertices of a tetrahedron, however, there are two ways to 
define its interior: either that part of S^ that includes the north pole (for 
example), or that part that doesn't. Taking one of these to be positive, say 
< y < 1 (in units of 27r2, the volume of S"^), the other will be y-1 < 0. Thus 
we cannot simply calculate the winding number of a given field configuration. 
We can resolve this ambiguity by defining the volume of a tetrahedron to 
satisfy always \V\ < 1/2. 

But now consider a tetrahedron whose volume is 1/2 — e. A small change 
in one of the vertices can push its volume to 1/2 + S, which will now be 
interpreted as —1/2 + S, while the volumes of its neighbors change hy 6 + e. 
Thus the winding number by our convention has jumped by 1: Topology has 
fallen through the lattice. The only way to avoid this is to force the field 
configurations to be sufhcicntly smooth that a tetrahedron volume stays well 
away from ±1/2. Wc can do this by choosing a sufhciently stiff kinetic term, 
such as Q 

^i = (a- l)^log(a„ -an+^i -a), (5) 

with [/ = flQ + icLiffi. This constrains a • a' > a. We find that setting a = 0.1 
keeps the tetrahedra well away from ambiguity. 

With the action (0) alone, solitons will collapse to radii on the order 
of the lattice spacing before they notice that the action isn't quadratic in 
derivatives (a la Derrick's Theorem) . To make sizable skyrmions possible, we 
add a Skyrme term 



<S'2=42J2J |(an+/i - an+p) (an+/i+p - a„) 
n fi>u 

- [(a„+/i " a„+j>) • (a„+/,+^ - a„)]^| . (6) 

The total action is S ^ (3iSi + (3282- 

Will the system still tunnel between topological sectors? In principle, 
yes. A non-local updating scheme could create a smooth skyrmion at a blow. 
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which would change the winding number by 1. This won't happen with a 
local algorithm, such as local Metropolis updating. 

This determines how to choose an updating scheme. Thus, if we are 
interested in doing thermodynamics with a chemical potential, summing over 
topological sectors according to 

Z(M)-^e^"Z„, (7) 

n 

we would choose a scheme that does nucleate smooth skyrmions (or even 
point-like defects). Our interest, however, is in studying the properties of a 
single skyrmion, so we trap it in our computer by insisting on local updating 
that preserves winding number. 

4 The Skyrme model in three dimensions, and the future 

As a first application, we have looked El at the three-dimensional Skyrme 
model, which may be thought of as a Ginzburg-Landau theory for the full 
quantum theory at nonzero temperature. The bare couplings in the action 
(at fixed cutoff) depend on physical parameters such as f^, the baryon mass, 
and the temperature, but we have not yet tackled this matching problem. We 
ran Monte Carlo simulations in the single-skyrmion sector and measured the 
mean-square radius {R'^) of the topological density defined above. We found, 
to our surprise, that our lattice action S admits a multitude of metastable 
minima, configurations of the skyrmion that have (generally) smaller radii 
than the skyrmion that actually minimizes the action. As the "temperature" 
of the 3d model — namely, l//3i for fixed /?2//3i — is raised, more of these lo- 
cal minima come into play and the skyrmion actually ends up shrinking as 
it is heated. Recall that in this Ginzburg-Landau theory, this temperature 
refiects both thermal and quantum fluctuations of the 4d quantum theory. It 
is possible that renormalization of the ratio /?2//?i will destroy this effect. 

Directions for further study include simulation of the full 4d theory and 
renormalization of its couplings (at fixed cutoff) ; study of the dependence of 
the size and shape of the equilibrium skyrmion on the physical temperature 
and density; and measuring the sensitivity of various results to the chosen 
form of the lattice action. 
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